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Abstract 

We study how to measure the current structure of the process that B 
meson decays into two unstable fermions fa and in model independent 
way. We use the momentum distributions of subsequent decay products 
affected by fafb spin correlation. We have found the following: (1) We 
can extract the absolute values of two effective coupling constants from 
the opening angle between the particles decayed from fa and /t,. (2) We 
can extract the real part of the interference from the energy distribution 
of one of the decayed particles from fa or /;,. (3) No new information can 
be obtained from the energy distribution of two decayed particles from 
fa and /(,. (4) The imaginary part of interference is extracted from the 
azimuthal angle asymmetry of final-state decay products. (5) If only one 
of two fermions is unstable, we can extract the real part of interference 
from each of the energy distribution and opening angle distribution. We 
show several simple examples. 
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1 Introduction 



A huge number of B mesons are produced in B-factories. They are used to 
confirm the Standard Model (SM). Almost of all the results suggest that the 
SM, and especially, Kobayashi-Mazkawa ansatz are reliable. Recently, we search 
for rare events and SM-forbidden phenomena in B-factories with high statistics. 
However, new physics has not been seen. 

To discover them, it is important to search through many modes and many 
physical quantities. They are, for instance, CP asymmetry, forward-backward 
asymmetry, left-right asymmetry, energy distribution, and angular distribution. 
We want to detect not only the decay width but also these quantities. Also, we 
want to analyze as many channels as possible using the unified form for sim- 
plicity, facility, and practicality. Another important thing to discover the new 
physics is making the reliable SM prediction especially for the non-perturbative 
QCD effect. Also for this purpose, determining many physical quantities is 
significant. 

In this paper, we consider the general B — ^ fafb decay modes, where fa and 
fb are arbitrary fermions and fa is the antiparticle of fa- The CP violation 
can be measured in some of these modes [1]. These modes can be divided in 
two types. One is the leptonic decay modes and another is the baryonic decay 
modes. The SM prediction in leptonic modes are [21 [3] 

Sr(B2 ^ r+T") ~ 2.8 x 10"^ 

Br(S° ^ r+r") ~ 8.9 X 10"^. 
The experimental upper bound is [1] 

Br{B'^ T+T-) < 4.1 X 10"^ (2) 

On the other hand, the branching ratios of baryonic modes are predicted as [S] 

Br{B° ^ S-A+) ^ 2.0 x 10"^ 
Br{Bi ^ S2A+) ^ 2.2 x IQ-^. 

The experimental upper bounds are for example, [3] 



Sr(S2 -> A" A) < 9.3 x 10"^ 
Sr(SS AA) < 3.2 x 10"'' 
Sr(SO^A-p) = (2.ll[j;7)xlO-5 

Br{B^ A-A+) < 6.2 x 10~^ (4) 



Br{B+ A+A) < 8.2 x 10"^ 

Br{B+ ^ SOA+)Sr(S^ ^ = (5.6i|J) x IQ-^ 

Br{B+ S°A+)Sr(S° ^ A'+tt") = (4.0 ± 1.6) x 10"^ 
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-6° — >■ Ajp and -B+ — )■ S°A+ have already seen. Also, some other modes 
are predicted to be seen in near future by the SM or other models. Comparing 
the experimental result with the model predictions of current structure, we 
try to discover new physics, select a reasonable model, and consider the non- 
perturbative QCD effects. 

The modes which decay into unstable particles decrease the efficiency since 
it is difficult to detect the events, however these modes have the advantage in 
correlation detection. The correlation is detected as momentum distribution of 
a and h, which are the decay products of fa and fb, respectively. When we deal 
with these modes, we have to consider the whole process of 

B ^fa + fb ^ b + anything 

I (5) 
' — > a + anything 

because we cannot detect the intermediate state fa + ft- 



1.1 By^fafb decay 

From the partially conserved axial current relation, the general Bq fafb decay 
amplitude is given by [2] 

A, = ^fBmnGF[{C!. + ^^^^^^^^C\){fbi,fa) + (C| + (/,/,)], 

TUB niB 

(6) 

where /s, rn_B, and Gp are B meson decay constant, B meson mass, and the 
Fermi constant, respectively; ma and mb are fa and fb masses, respectively; C|>, 
C|, C^, and Cy are the complex coefficients of pseudo scalar, scalar, axial, and 
vector currents, respectively; The superscript q represents the valence u, d, s or 
c quark in B meson. 

In charged B meson decays, we simply set 

(7) 

niB 

On the other hand, in neutral B meson decays, considering the i?*^ — B'^ mixing 
effect, we set [5] -[5] 

|i3°(i))=.9+(i)|S") + VW|S°), 

(8) 



I _|_ g-iAmstgiArBt 



where; t is the time started when 5° is created; q/p is the ratio of 5° to 5° 
in B° mass eigenstate; is the B^ total decay width; Atub and AF^ are 
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the mass deference and decay width difference between heavier and lighter 
mesons; Hence, the time dependent effective amplitude takes the form 



where 



Ag{t) = ifBrnsGF [Ci {fbjsfa) + CMa)] , (9) 
5i = {5+(i)Ci + V(i)<5i}, 



C2 = {g+{t)C2 + ^9-{t)C2}. 



(10) 



These parameters appear in the differential decay width in the form of |Ci|^, 
IC2P, -Re[CiC|], and /m[CiC|]. These quantities depend on the decay time. 
However, for \p/q\ = \Ci/Ci\ — \C2/C2\ = 1 and ATb = 0, integrating over 
the time and summing over decays and decays, these quantities be- 
come_s 2|Ci|VrB, 2|C72|VrB, iRe[C\C^] + Re[CiC^]) /T b , and {Im[CiC^] + 
Jm[CiC|])/rB, respectively. Hence, we omit the time dependence and the 
tildes on Ci and C2 in most of the rest of this paper. 

We want to give the B faft partial decay width in which fa and have 
particular polarizations. Thus, we introduce the polarization vectors and 
of fa and fb, respectively. These vectors have the constraints (s")^ = (s'')^ = — 1 
and s^-kf^ = s^-kf^ = 0, where fc/^ and kf^ are /„ and /& momenta, respectively. 

In B rest frame, the differential decay width oiB ^ fa{s°')fb{s'') is given by 



dT flGp. 



1 1 In. 1 1 oh ^ 



(11) 



where 



IPI 



^J (m\ - {nig - m.bY){m\ - {ma + mbY) 
2mB ' 



|2 m| - (m„ -mb)^ |2 "^| - (m„ + m6)^ 
Ui = \Ci\ I-IC2I , 

m% - {mg - TUbY , 1^ |2^7i| - (wg + mb)^ (12) 



D2 = -|gi r"^ + 1C21 

D4 = -2Re[C\C;]mamb-fg-fb{l3a + Pb), 
D5 = -2Im[CiC2]m,amb^alb{Pa + M, 



1,0 ' 



la 



lb 



1 



1 

V^b 



ma 

h 

nib ' 



(13) 
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and n is the solid angle of kf^. 

The general B ^ fafb ^ a + 6 + anything differential decay width is written 

as 

dT ^ ^ dr{B ^ /a(s°)/fc(s^)) 

dn d^ka d^kb ^ dt dn 

is", is' 

dBr{fa{s'') a + anything) (24) 

¥ka 

dBr{fb{s^) -^h + anything) 
dFkb ' 

where S impUes sum over polarizations, ka and fcf, are the momenta of the 
particle a and h in fa and fb rest frame, respectively. The differential branching 
ratios of fa and fb are written in Appendix [A] 

In writing the explicit form of the decay width, we will use the following 
notation (Sec Fig. [TJ): In B rest frame, fa is oriented in the positive z-axis 
direction. The zenith angles of a and h directions in B rest frame are Oa and 
9b^ respectively. The azimuthal angle between a and b directions is (p. dz is the 
distance between fa and fb decay points. 




Figure 1: B meson rest frame, fa is oriented in the positive z-axis direction. The 
zenith angles of a and b directions are 6a and 6b, respectively. The azimuthal 
angle between a and b is 0. The distance between fa and fb decay points is dz. 
a and b have the momenta and k;,, respectively. 
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In the massless limit of a and 6, we obtain the general formula 
dV 



where 



dyadfla dybdflb 

2 2 f2 r<2 

L ma mt J 

— {D2 sin 9a sin 6}, cos (j) — Di cos 9a cos 6'b + D5 sin fl^ sin 9^, sin (/>) 6*2(1/0) G2 (j/b 

(15) 

i?ra = Br{fa -> a + anything), 

Br6 = Br(/5 ^ 6 + anything), ^^^^ 
2Ea 2Ei, 

Va = , yb = — , 

TUa rUb 

where Ea and E^ are a and b energy in fa and ft rest frames, respectively; 
G° 2(2/0, fc) are the functions which are defined in Appendix El The massless 
condition of particles a and b are reasonable because most of r decay into /i, 
e, or pions, and substantial unstable baryons decay into a lighter baryon and 
pions, photons, and/or leptons. They have at most about 100 MeV masses, 
which are enough smaller than the masses of r and any baryons. 

Using the general formula (fT5|). we first derive the partial decay width. In- 
tegrating over dyadftadybdfli,, we have 

T = BraBrb^\p\D,. (17) 

This width contains the factor Di. We determine this coefficient, first. How- 
ever, we want to know the relation between \Ci\ and IC2I. Moreover, we want 
to know how is the relative phase between Ci and C2. That is what we will do 
in this paper. 

This paper is organized as follows: In Section 2, we consider the energy dis- 
tribution of a to determine Re[CiC2]- In Section 3, we consider the distribution 
of opening angle between a and b to determine \Ci\ and IC2I, separately. In 
Section 4, we consider the azimuthal angle asymmetry of a and b to determine 
Im[CiC2]- In Section 5, we discuss the case that fb is a stable fermion. In Sec- 
tion 6, we show some examples of baryonic mode. In Section 7, we summarize 
our analysis. 



2 Energy Distribution 

In this section, we study the energy distribution of the final-state particle a or b. 
For definiteness, let's say that we want to investigate the a energy distribution. 
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The prescription to derive the energy distribution formula in B rest frame is 
as follows [9]: First, we multiply the delta function S[xa — 2/a(l + /3a cos0a)/2) 
by Eq. (ITS]) , where Xa ~ E'^/Ef^, and Ef^ and E'^ are fa and a energy in B 
rest frame, respectively. Xa means a normalyzed energy of particle a in _B rest 
frame. Next, we integrate over dyadfladybdflb- Then, we have 



Tl^= I dy^^{yaG'i{ya) + ^^{2Xa-ya)G'^2{ya)]. 
T dXa J PaK ^ Di maPa i 

Here, J dya means 



(18) 



dya = j dya0[xa]6\ } - a;Q] + j dyaO[xa - ^—^]0[^-^^^ - Xa]- 

TTWZ TT^ ^^^^ 

The expression (jl8p suggests that a energy dependence can be used to determine 
the coefhcient D4, which contains _Re[CiC|]. We note that b energy dependence 
can also be used to determine D^, similarly. However, no new information is 
obtained by the energy distributions of both of a and 6, namely, dV / [dxadxb). 
This is because D2 and D5 terms in the general formula (|15p vanish when we 
integrate over the azimutlial angle 

2.1 Example 1 - r+ Decays into /i+z/^P^ 

As a simple example, we calculate /i"*" energy distribution of 

^ +^ 4.- 
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In this case, we can set G^{ya) = 3 - 2j/p, G^{ya) = 2y^ - 1, Xa = Xb = 5, and 
Pa — Pb — — Am1/m?^ = /3. Hence, we have 

1 dr 

r dxf^ 



813x1 (9(l-/32)_4(3 + ^2)^^) 
3(1-/32)3 

IC1I2 + IC2P/32 3(1-/32)3 J ^^fJ 2 



+ 



• 5(1 + /3)3- 4(9(1 + /?)- 8x^)0;; 



+ 



6(l + /3)3 

2j?e[CiC2l ((1 + Pf ~ 12(1 + 2/3)(l + /3)2:2 + i6(l + 3l3)xl) 
|5i|2 + |C2|2/32 6(l+/?)3 



(21) 



Here, we put tildes on Ci and C2. Integrating F and dT/dXf^ over the time, 
and summing them over B'^ decays and 5° decays, the energy distribution is 
represented by Eq. ^ replacing |Ci|2, jCap, and Re[CiC^] with |Ci|2, jCzp, 
and {Re[CiC^] + Re[CiC^])/2, respectively. 

We depict this distribution and perform a Monte Carlo simulation (MC) to 
estimate the error of {Re[CiC^] + Re[CiC^])/\2CiC2\ in Figs. [H 

Fig. [5] represents the |Ci| = IC2I = 1 case. Similarly, Fig. |3]and Fig. H] rep- 
resent the {|Ci|,|C2|} = {1,0.1} and {|Ci|,|C2|} — {0.1,1} cases, respectively. 
The interference effect emerges when \Ci \ ~ |C2|/3. We note here that in these 
figures, Re[CiC2] independent point where x^ ~ 0.4 is one of the solutions of 
identity, (1 + /3)3 - 12(1 + 2/3) (1 + l3)xl + 16(1 + 3l3)xl = 0. 

In Fig. m the MC is performed in a sample of 2000 events for {Re[CiC2] + 
Re[CiC2]) /\2CiC2\ — 0. The number of events are given as follows: The Super 
KEKB will make about 50 ab~^ integrated luminosity. The e"'" + e^ — > T(45) — >■ 
^d^d cross section is about 10~33 cm2. The — r+r~ branching ratio is in 
Eq. ID). The r+ ^^v^Vt branching ratio and the r+ — > c^v^Dt branching 
ratio are about 0.174 and 0.178. These are essentially the same events for 
the massless limit of daughter fermions. Therefore, about 2000 events will be 
available. The efficiency of this mode is in fact very low. However, we here just 
ignore it. The MC result is {Re[CiC^] + Re[CiC^]) /\2CiC2\ = 0.11 ± 0.10. 

In Figs. [3] and [H the MC is performed in a sample of 20000 events for 
{Re[CiC-i]+Re[CiC^])/\2CiC2\ = 0. The MC results in {Re[CiC^]+Re[CiC^])/ 
I2C1C2I are -0.05 ±0.21 and 0.17 ±0.12, respectively. 

Jn the SM, |Ci| > IC2I [2]. Then, the coefficient 2Re[CiC^]/ {\Ci\'^ + 
/3|C'2p) is nearly zero. This situation is also realized if we set Re[CiC2] = 
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2000 events 

dr 
rdx, 

3.0 r 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: The /i energy distribution of B'^ — > r+r^ and subsequently — >■ 
/u"'"(e"'") + f^(e) + or t~ — >■ ^~{e~) + Vfj_(f.)+VT decay, and their CP conjugate. 
The horizontal axis is the normalized ^ energy a;^. The vertical axis is the 
time integrated differential decay width dT /dx^ over the time integrated partial 
width r. We set [Ci | = IC2I = 1. The solid line, dashed line, and dot-dashed line 
represent {Re[CiC^] + Re[CiCl\) /\2CiC2\ = {0,1,-1} case, respectively. The 
MC result in a sample_of 2000 events for {Re[CiCi] + Re[CiC-^]) /\2CiC2\ = 
is {Re[CiC^] + Re[CiC^])/\2CiC2\ = 0.11 ± 0.10. 



Re[CiC2] — 0. Therefor, the shape of distribution is the same as {Re[CiC2] + 
Re[CiCl)/\2CiC2\ = case. 

B — >■ T^T~ energy distribution is very interesting since we can investigate 
the current structure of new physics. 

For instance, Ref. 10^ expresses the Higgs induced operators for the tran- 
sition b — )■ sfi^ fjT . It is easy to transform them for the transition h — !■ dr^r" . 
Concletely, it is realized by the deformations — > rrir , F23 F13 , ^"3*2 — > F^i , 
^j, ^ T, and s — > d. After that, the coefficients Ci and C2 in this paper are given 

by 



1 Tot- sin /? 

(22) 



^'--2G^2vcosl3 Ml^^'' + ^^' 



1 Tot I sm{a — /3) cos a cos(a — /3) sina , , ^, ^ 

= ~2G^2^ [ Ml Ml ' - 



In Eq. ([22|) . v, 13, a, Mh, M^, Ma, F13, and F31 are vaccum expectation value, 
Higgs mixing angles, hevier nutral Higgs mass, lighter nutral Higgs mass, CP 
odd Higgs mass, and b — d coupling constants, respectively, which are defined 
in Ref. [TU] . This contribution can compete with or even dominate the SM one. 
Especially, in Ma — > 00 limit, Ci contains only the SM effect and C2 contains 
only the 2HDM contribution. 

On the other hand, the supersymmetric SM (SUSY) models without R-parity 
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20 000 events 

dr 
rdx, 

2.5 r 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3: The absolute values of coefficients are \Ci\ — 1 and IC2I = 0.1. The 
number of event is 20000. Others are the same as Fig. [2l The MC result is 
{Re[CiC^]+ Re[CiC^])/\2CiC2\ = -0.05 ±0.21. 



[TT] suggest the coefficients 



C2 



Ik 



2tot A'Jj.^Al 



3k3 3kq 



[CI] 



Ik 



(23) 



{k=£3} 



where Xijk and A'^j. are the coefficients of LiLjlf. and LiQjd'f. type couplings, 

respectively, where L, l'^, Q, and df^ are the lepton doublet, lepton singlet, quark 
doublet, down-type quark singlet superfields, respectively; m,-. and rriq^ are slep- 
ton and squark masses, respectively. 

Moreover, leptquark models [T^], the topcolor-assisted technicolor model 
[13], and the Babu-Kolda model [Tl] also deform the SM energy distribution, 
while the energy distribution in the multiscale walking technicolor model [15] 
is the same as the SM one. Considering the ratio Br{B TT)/Br{B — > /i/i), 
SUSY models without R-parity [TT], leptquark models [T^j, and the topcolor- 
assisted technicolor model [13] predictions differ from the SM one, while the 
Ref. [TD], the Babu-Kolda model [T3], and the multiscale walking technicolor 
model [ini predict the same value as the SM one. These charasteristic features 
of models are available to distinguish them. 

All of these models predict that the Br{B tt) may become larger than 
the SM one. 

2.2 Example 2 - t+ Decays into -k^Vj. 

We show here another example, in which B^ decays into r+r^ and subsequently 
r+ ^ 7r+ + Dr. In this case, we can set Gl{ya) = S{1 - yn)/yl, G2{ya) = 
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20 000 events 

dr 
rdx. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 4: The absolute values of coefficients are \Ci\ — 0.1 and IC2I = 1. The 
number of event is 20000. Others are the same as Fig. [21 The MC result is 
{Re[CiC^]+Re[CiC^)/\2CiC2\ = 0.17 ±0.12. 



-6{l - yT,)/yT,, \a 1, and Pb = ^1 - 4m2/m| = /3. Hence, the energy 
distribution is 



In Figs. [S][71 we depict the time integrated distributions and the MC re- 
sults for {Re[CiC^] + Re[CiC^]) /\2CiC2\ = as the previous example. Fig. 
[5] represents the \Ci\ — IC2I — 1 case. Similarly, Fig. |6] and Fig. [7] represent 
the {|Ci|,|C2|} = {1,0.1} and {|Ci|,|C2|} = {0.1,1} cases, respectively. This 
mode is more suitable to understand the — ?► r+r" current structure than 
preceding one, since two-body decay does not dilute the polarization unlike the 
previous case, even though the t'^ tt"*" -|- Vt branching ratio is about 0.11, 
which is smaller than the previous case. 

The results of MC are as follows: 

Re[CiC^] + Re[CiC^] 
2IC1C2I 

-0.11 ± 0.07 for \Ci \ = IC2I 1, 600 events (25) 

-0.21 ±0.15 for jCil = 1, IC2I = 0.1, 6000 events 

-0.036±0.085 for |Cil = 0.1, 1(721 = 1, 6000 events 



3 Opening Angle Distribution 

Here, we consider the opening angle Q between particles a and h 'm B rest 
frame. The prescription is similar as Section [2j however, this time we multiply 
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600 events 




Figure 5: The 7r+ energy distribution of B° — >• r+r" and subsequently r"*" — >■ 
77+ + Pr or T~ Tr~ + fr decay, and their CP conjugate. The horizontal 
axis is the normahzed 7r+ energy a;^. The vertical axis is the time integrated 
differential decay width dT/dx.^ over the time integrated partial width F. We 
set |_Ci] = IC2I = 1. The number of events to perform MC for {Re[CiC^] + 
Re[CiC^])/\2CiC2\ = is 600._ The solid line, dashed line, and dot-dashed 
line represent {Re[CiC2] + Re[CiC2])/\2CiC2\ = {0,1,-1} case, respectively. 
{Re[CiC^] + Re[CiC^])/\2CiC2 \ = -0.11 ± 0.07. 



the different delta function 

sin 9a sin 9b cos (j) - ialb{Pa + cos 9a) {Pb - cos 9b) 



cos© — 

Then, the result is 
1 dV 



(1 + /3a cos 6'a)(l - Pb cos 6b)^alb 



FdcosO 

1 



BMax(e,e„) 
d cos 9a I d cos 9b 



[1+ fia C0S^a)(l - fib COS 9b)-falb 



(26) 



7a76 {(1 + /3a COS 9a){l - f% COS 9b) COS 6 + (/3a + COS 9a)il3b - COS 9b)y 



sm 



, sm 



|l + ^(— COS 9a{G^2{ya)} + — COS 9b{GUyb))) 

I. JJi rua nib 



D2 

[Yrlo.lb{{'^ + PaC0s9a){l - /36 COS ^fe) COS 6 + (/3a + COs6la)(^6 - C0s6b)} 



-cos9aCos9b\{Gl{ya))^G\{yb))\. 



(27) 
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6000 events 




Figure 6: The absolute values of coefficients are \Ci\ — 1 and IC2I — 0.1. 
The number of event is 6000. Others are the same as Fig. [S] {Re[CiC2] + 
Re[CiC^])/\2CiC2\ = -0.21 ±0.15. 



where we set 



(28) 



B M ax ,min 

_7a7b {(1 + Pa COS ea)(3b COS 9 + (/gg + COS 0a)}{{l + pg COS Og) COS 9 + (^gg + COS ggj/^b} 

7a76{(l + Pa cos9g)Pb COS 9 + {/Sg + cosea)}^ + sin^ Og 
sin 6lg ^ sin^ e^g + 7a{(/3a + cos 0^)2 - (1 + /3a cos 6'a)2 cos2 9} 



± 



7a7b {(1 + /3a cos ea)Pb COS 9 + + COS 6lg)}2 + sin^ Og 



(29) 



This expression suggests that the opening angle distribution determines \Ci\ 
and I C2 1 , separately, via the coefficients D2 ■ 

If fa = fb, and the decay modes of fa and /& are the same, for example, 
^ T+r^ — !■ TT+TT"/?,-;/,- mode, the second term in Eq. (P7| which has 
the coefficient D^/Di will vanish because this term is antisymmetric about 
the cos 9a + cos 9b = line, on the other hand, the domain of integration is 
symmetric. 

3.1 Example 3 - t"^ Decay into fj."^ 

We here show a simple example that B'^ decays into t^t^, and subsequently, 
they decay into /i+ + /i^ + + i?^ + + Vr- In this case, we set (62(2/0)) = 



12 



dr 

rdx. 



6000 events 



1.5 - 




1.0 - 



0.5 - 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 7: The absolute values of coefficients are \Ci\ — 0.1 and IC2I = 1. 
The number of event is 6000. Others are the same as Fig. [S] {Re[CiC2] + 
Re[CiC^])/\2CiC2\ = -0.036 ±0.085. 

{G\{yb)) — 1/3. The numerical result is depicted in Fig. [8] We perform the MC 
for (|Ci|2 - |C2p/?2)/(|Ci|2 + = in a sample of 35000 events, which 

corresponds to the 100 times of 50 ab'^. The resuh is (|Ci p - |C2p,32)/(|Ci 1^ + 
IC2P/32) = -0.15 ±0.18. 

In this figure, the increase near cosG = — 1 is caused by the back-to-back 
Lorentz boost of fa and /{, along the z axis. 

If new physics affects this mode substantially, we may detect the distribution. 
In that case, B — > r"'"r~ opening angle distribution is usefull to distinguish new 
physics models. In the SM, \Ci\ » IC2I [2]- Then, the shape of distribution 
is the same as (|Cip - |C2p/3^)/(|Cip + |C2p/3^) = +1 case. However, many 
new physics models, for example, Refs. [TU], [H], [12] j [13]) and [Hj, deform the 
shape of distribution near cos 9 = — 1. 

3.2 Example A - Decay into tt^ 

We show another example that 5° decays into r+r", and subsequently, they 
decay into 7r+ + 7r~ ^ v^- -\- i>T- In this case, we set (Gf (j/a)) — {G\{yb)) = 1. 
The numerical result is depicted in Fig. [Qj We perform the MC for (|Ci|^ — 
|C'2p/3^)/(|Cip + |C2p;52) = in a sample of 40 events, which will given in 
the Super B-factory The MC resuh is (jCip - \C2\^ p^) / (\Ci\^ + \C2?P^) = 
0.03 ±0.49. 

In this figure, the dashed line decreases near cos 9 = — 1, on the other hand, 
the dot-dashed line increases there. The reason is as follows: In this case, the 
second term in Eq. J^Tj) is vanished, and we can set fia — Pb = P- Therefore, 
for cos 9 = — 1, Eq. (P7|) is proportional to 




(30) 
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35 000 events 
dr 



rdcos© 




Figure 8: The opening angle distribution between /i'^ and fi" in B'^ t'^t^ — >■ 
11^ + + VT-v-yVfiVyi, decay. The horizontal axis is cos 0, and the vertical axis is 
time integrated dF/dcosO over time integrated F. The solid line, dashed line, 
and dot-dashed line represent {\Cxf - |C2|2/32)/(|Ci|2 + \C2\^P'^) = {0, 1, -1} 
cases, respectively. The MC result for (|Ci|2 - \C2\^ 0^) j {\Cx\^ -h jCzp/J^) = is 
(|Ci|2 - |C2p/32)/(|Ci|2 = -0.15 ±0.18 in a sample of 35000 events. 

The factor 1 -|- D2I D\ becomes zero when C2 is zero, and it becomes maximum 
when C\ is zero. 

4 Azimuthal Angle Asymmetry 

Generally, the trajectories of a and h draw the skew lines since ja and /b have the 
finite lifetimes. If the vertex detector of B-factory could detect the decay points 
of /a and /fc, we were able to determine (/) dependence of dF, and then Im^xC^. 
However, some of /a and decay into one-prong modes, the polarization effect 
is diluted in the many body decays, and/or the vertex detector does not have 
sufficient resolution to detect the decay points accurately. Thus, we consider 
another method to determine I'm\C\C'2\- 

Since is the azimuthal angle between a and h as depicted in Fig. [1] the 
Lorentz boost along z direction has no effect on this angle. Thus, the delta 
function is unnecessary unlike the Sections [5] and [H 

The trajectories of a and 6 in B rest frame are written as 

qa(ia) = {sin 6*0^0,0, (/3 

qb(^b) = {sin6'bcos0i6,sin6'bsin0i6, (-/3b7b -|-7;,cos6'f,)4}, 

where ta and are the parameters. 

The vector product of qa(io) and qbih) for dz 0, ta > 0, and tb > takes 
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40 events 
cir 

rdcos© 
1.5 - 




cos© 



Figure 9: The opening angle distribution between tt"^ and vr^ in i?*^ — r+r^ — >■ 
TT"*" + 7r~ + VtVt decay. The horizontal axis is cos0, and the vertical axis 
is dV /iVdcosQ). The solid line, dashed line, and dot-dashed line represent 
(|Ci|2 - |C2p/3^)/(|Ci|2 + |C2|2/32) = {0,1,-1} cases, respectively. The MC 
result for (jCip - |C2|2/32)/(|Ci|2 + \C2?P^) = is - |C2|2/32)/(|Ci|2 + 

|C2p/3^) = 0.03 ± 0.49 in a sample of 40 events. 



the form 



(~la{l3a + COS 6'a) sin 6*6 sin0 
la{Pa + COS 6*0) sin 6*6 COS0 - SU\9a^b{-Pb + COS 6*6) 
sin 9a sin sin (f) 




Meanwhile, the difference between (ia{t'a) ^-^d qb(t[,) is 

(sin Oat'a — sin ^j, cos (fyt'f^ \ 
-sin 0b sin , (33) 

la{Pa + COS 61^)^^ +d,- -fbi-Pb + COS 61,,)^;, / 

where t'^ and tj, are the parameters of (la{t'a) and cib{t[,)- tj, and tj, take arbitrary 
values. The scalar product between Eq. p2p and Eq. (p3|) is given by 

qa{ta) X qh(tfc)|d,-i-o,t„>o,tb>o ■ {qa(il) - <ibit'b)} = tatfodz siu 6*0 sin 6*6 s'mcf). 

(34) 

This quantity becomes plus as < (/) < tt and minus as tt < </> < 27r. The sign 
of sincj) is determined event-by-event (See Fig. [TU)) . Then, the azimuthal angle 
asymmetry 

TT 
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Figure 10: Interpolating the trajectories, the skew hnes take two types of ahgn- 
ments. The left one corresponds to < cj) < tt. The right one corresponds to 
TT < < 27r. 

gives us the coefficient D5, which is proportional to Im[CiC2]- 

We perform the MC for B° — > r++T^ and then t+ — > tt+P,- and ^> tt^Vt-- 
We set (/m[CiC2]+/m[C'iC'|])/|2CiC2| = and summarize the results in table 
1. 



Table 1: The MC for {Im[CiCI]+ Im[CiCi]) /\2CiC2\ = and some conditions 



\Ci\, IC2I 


number of events 


Irn[C^Ci\ + Irn[C^Ci\ , 


= IC2I - 1 


40 


-0.40 ±0.64 


= IC2I = 1 


400 


-0.053 ±0.19 


|Ci| = 1, IC2I =0.1 


4000 


0.078 ±0.39 



5 Only One of Two Fermions is Unstable 

If fb is a stable particle, for example, B'^ -)■ A~p, 5+ Sc(2455)''p, 
t+i/t-, S° — >■ T+e~, and the general formula ([T5]) is modified to 

form 

t4o- = ^^a/^4^|p||i?iG?(j/„) +i?4— cos0,G^(y,)|. (36) 
dyadila 47rAa 27r L J 

Then, by the similar calculations, the partial decay width is 

r = i?r, 4^|p|A, (37) 

the a energy distribution is 

f dya--^{yaGUya) + ^^i2Xa-ya)G^2(ya)}, (38) 
1 dXa J PaK ^ L>1 rUaPa J 
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where / dya means the same as before, and the distribution of the opening angle 
0' between and a is 

1 dr _ 1 1 - /3q r DiTUb cos 9' + Pa 



ll _ ^4mfc cosu +Pa ^ 

rdcos9' 2 (l + /?aCos9')2 1 i?imal + /3aC0s9'^^^^"^^i' ^ ' 

where 

, Pa + cose a 
cos 9 = — . (40) 

l+Z^aCOS^a ^ ' 

Both of these two distributions give D^. These are used for a cross-check. 
However, we cannot pull out D2 and D^. The energy distribution is useful even 
if fh is a missing fermion except for neutrinos. If fb is a neutrino, rub — > and 
the second terms in both of Eqs. ((38|) and ((39|) are vanish, and then we cannot 
determine 1)4. 



5.1 Example 5 - S° r+/x ,r+ tt+z/^ 

We consider the lepton flavor violating — > T+/Lt~ decay, and subsequently 
r+ decays into tt+Pt-. In this case, we can set Gi{ya) = S{1 — y)/y^, ^2(2/0) = 
-d{l-y)/y, Xa^l, and (G^(ya)) = 1. 
The a energy distribution is 



1 dr _ 1 r 

r dXa Pr ^ 



l-^^(2x.-l)}%.-^^M^^-x.]. (41) 



Di rrirPT 

On the other hand, the opening angle distribution takes the form 

1 dV I l-pl ( Diui^ cos9' + /3. ^ 
rdcos9' 2 (l + ;3^cos9')2 I Di 1 + cos9' /' ^ ^ 

These distributions and the MC for {Re[CiC^] + Re[CiC^])/\2CiC2\ = in 
a sample of 10000 events are depicted in Figs. [TT] and [T^l The results are 
{Re[CiC^] + Re[CiC^])/\2CiC2\ = 0.27 ± 0.29 and -0.03 ± 0.30, respectively. 
They have almost the same errors. 

The SUSY models without R-parity ^16| suggest the coefficients 



4G 



(43) 



where C^g, C^, and Cy are defined in Ref. 16 as 

cts-Y.'^' ct^^Y.'^' ^i-Y.^^. («) 
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10 000 events 




0.2 0.4 06 0.8 1.0 



Figure 11; The energy distribution of 7r+ in B*^ — ^ T+/i , t+ — s> tt^Vt for 
It'll = |C'2| = 1- The sohd line, the dashed hne, and the dot-dashed hne 
represent (Re[CiCi] + Re[CiC^])/\2CiC2\ = {0,1,-1}, respectively. The MC 
is performed for {Re[CiCi] + Re[CiCj])/\2CiC2\ = in a sample of 10000 
events to resuh in {Re[CiC^] + Re[CiC^])/\2CiC2\ = 0.27 ± 0.29. 



where rrij}. is squark mass. In this model, the branching ratio can become about 
10"^, which is the same order as the experimental upper bound, and it has a 
relation, Br{B — >■ /i/x) < Br{B — > r/i) < Br{B -> rr). Moreover, since this 
mode has only one neutrino, the efficiency is much higher than that of i? tt 
mode. 



6 Example 6 - 5+ ^ E^K^ and A+p Decay 

Now we show two baryon modes for examples. The first one is 5+ '^c-^t 
decay. According to Ref. [5], in the SM, we have the relation for £?+ 

= \2 = 0.10, (45) 

where m^^ and are S[? and masses, respectively. Then, we predict 

^=0.89, ^ = -0.45, ^^0. (46) 

Di ' Di ' £>! ^ ' 

This prediction is available for the test of Ref. [S] . 

The next one is B^ A+p decay. According to the Ref. [17], which uses 
the factorization, 

^=0.34. (47) 
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10000 events 
dr 




Figure 12: The distribution of opening angle between /i and 7r+ in the same 
process as Fig. [TTJ The solid line, the dashed line, and the dot-dashed line 
represent {Re[CiCI] + Re[CiC^])/\2CiC2\ = {0,1,-1}, respectively. The MC 
is performed for {Re[CiCi] + Re[CiCj])/\2CiC2\ = in a sample of 10000 
events to resuh in {Re[CiC^] + Re[CiC^])/\2CiC2\ = -0.03 ± 0.30. 

Then, if C1C2 has no relative phase (or i?e[CiC|] — IC1C2I), 

|i = -0.52. (48) 
On the other hand, according to Ref. [18], which uses the pole model, 

^^0,77 (49) 

in B''^ —i' A+p decay. Therefore, if C1C2 has no relative phase (or Re[CiC2] = 
IC1C2I), 

|i = -0.90. (50) 

These two models suggest different branching ratio and IC2/C1I. The dif- 
ferent branching ratio may be corrected by the non-perturbative QCD effect. 
However, IC2/C1I eminently represents the feature of each model. Hence, we 
can test which model works better. 

The QCD effect may pollute new physics effect. However, at least, one of 
these observables is measured to considerably differ from (|46l) . or (|15)) and ([SO]) , 
we should take into account new physics. We point out that more observables 
are desirable for new physics discovery. 

Theoretically, deriving the precise expressions of Gi and G2 is not easy. 
However, we are interested in i? — > faft decays. We assume that there is no 
new physics contribution in the and decays. So, it is not necessary for 
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us to know the expressions of Gi and G2 theoretically if we can determine it 
from the experimental data. Actually, Refs. [4], [19], [20], and [21] suggest the 
Ac polarization, also Ref. [55] measures the E!° polarization. In this reference, 
S° decays into S~7r+. The two body decay makes Gi and G2 trivially. Then, 
these decays are given by 



TV ^ 

dBr{A+ ^ 7r+ + A) _ 2Sr^+ 



■ — s ■ k^- 



2 ^"AoS " ■ K+ 



(51) 



where Br^o = Br{E^ — > n~ and Br^+ = -Br(A+ 7r+ + A); m^± are 

the charged pion masses, y^ri are normalized charged pion energies defined by 
the same manner as in Eq. (jl6p : and a^^ are the decay parameters; s"<= 
and s^<: are the polarization vectors of and A+, respectively; k^± are the 
TT^ momenta, respectively. 

Fig, [13] explains the B+ + A+, A+ A + 7r+, and S° p+ + 

7r~ decay Monte Carlo simulation with 320 and 16000 events. The horizontal 
region is determined by the tt'^ energy distribution and the diagonal region 
is determined by the 7r"'"7r~ opening angle distribution. The dot represents 
Eq. (|46p . The simulation is performed with this parameter set and the decay 
parameters, as.^ = —0.6 and = —0.91. We have to consider that as^ has a 
large ambiguity as.^ = —0.6 ± 0.4. So, we explain as.^ dependence in Figs. |T3] 
andHS] 

Fig. [Tn]represents the 7r+ energy distribution of B^ — )■ A+p and A+ — )■ A+tt"*" 
decay chain. The hghter gray explains |C2|/|Ci| = 0.77 case and the darker 
gray does |C2|/|Ci| = 0.34 case. These distributions have breadth caused by 
the ambiguity in = —0.91 ± 0.15. The results of Monte Carlo simulation 
with 200 and 10000 events are 

"Ac ^ -0.72 ±0.27, 
"Ac TT =0-83 ±0.04, 

respectively. 



7 Summary and Discussion 

We studied the current structure of S — >■ fafb decay modes using polarization 
effects. This can be applied to both of leptonic and baryonic decays, also, to 
the charged and neutral B mesons. 

The a energy distribution gives Re[CiC2]- If we consider no or small relative 
phase between Ci and C| , we can estimate the ratio of | Ci | and | C2 1 . The energy 
distribution of a and b gives no more information. 
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320 events 



16 000 events 




Figure 13: The allowed region determined by the energy distribution and open- 
ing angle distribution of -B+ ^ S° + A+ , A+ — > A + 7r+ , and S° ^> p+ + vr^ 
decay chain. The left and right figures are the results of MC simulation with 
320 and 16000 events, respectively. In each figure, the dot means the SM pre- 
diction. The horizontal lines explain the allowed region which is determined by 
the energy distribution. The diagonal region is allowed by the opening angle 
distribution. = -0.44 ±0.10 for 320 events. Di/Di = -0.44 ±0.01 for 

16000 events. 




Figure 14: Same as Fig. [13] but, we here use a^^ = —0.6 — 0.4. The diagonal 
region in Fig. [T3] becomes narrower and more vertically. 
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320 events 



16000events 




Figure 15: Same as Fig. [13] but, we here use as^ = —0.6 + 0.4. The diagonal 
region in Fig. [T3] becomes wider and is inchned more horizontally. 



200 events 10 000 events 

dr dr 

rdx, rdx, 

2.5 r , r 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Figure 16: The 7r+ energy distribution of K'^p and A+ — > A + 7r+ 

decay chain. The lighter gray explains IC2/C1I = 0.77 case and the darker 
gray does IC2/C1I = 0.34 case. These have breadth caused by the ambiguity in 
= 0.91 ± 0.15. The dots with error bar are the result of MC simulation for 
IC2/C1I = 0.77 case with 200 (left figure) and 10000 (right figure) events. These 
suggest that aA^DjDi = 0.72 ±0.27 (left figure) and aj^^D^/Di = 0.83 ±0.04 
(right figure). 
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The opening angle Q distribution gives \Ci\ and IC2I, separately. With the 
energy distribution, this gives us the relative phase between Ci and C2 up to a 
binary ambiguity. 

The azimuthal angle (j) asymmetry gives Im[CiC2]- We cannot detect the 
decay point in the one-prong events. Then, we cannot determine the cj) distri- 
bution. However, we can determine that 4> is larger or smaller than tt. This is 
enough to give Im[CiC2]- 

If one of two fermions is stable particle, we cannot determine D2 and D^. 
However, D4 is determined by each of the cos 8' and the Xa distribution. 

We predicted D4/D1 and D2/D1 of the baryon modes 5+ — > ScA+ and 
B° — >■ A+p. They are summarized in Figs. [T^lfTOl 

In the Examples 1-4, we derived the — s> r+r^ sample number ignoring 
the efficiency. Here, we try to consider it. According to Ref. [23], they conclude 
Br(B2 T+T-) < 4.1 X 10-3 using (232 ± 3) x 10^ data sample, which cor- 
responds to 210 fb"^. Hence, we need 7.9 x 10^^ (7.2 x lO'^ ab^^) data sample 
to discover a — !• r+r" event with the same efficiency as Ref. [53]. This 
efficiency can be improved, for example, by the semileptonic tagging method 
[24] . However, it is difficult to detect this mode in the SM case. If — > r+r^ 
is detected, it must be induced by new physics. Then, our analysis is useful to 
determine its current structure. 

In the neutral B decays, Ci and C2 are the functions of t as defined in Eq. 
PH)) . If we determine |Cip, IC2P, i?e[C'iC|], and Im[CiC2], respectively, then. 



using their t dependence, we can derive the time independent coefficients, 
result is summarized as follows: 



The 



IC2P ^ IC2I 



|C'ip,i?e[Ci^Ci*],/m[Ci^(7*] 
p* p* 

\C2\^Re[C2^C*2],Im[C2fy*2] 
p* p* 



Re[Ci — C;] + Re[-CiC;],Im[Ci—C;] - Iml^CiC;] 
p* p P* P 



Re[Ci^C;] - Re[^CiC;],Im[Ci^C;]+Im[^CiC;] 



(53) 



To determine these quantities, we need a huge number of statistics. 

In Examples 1-5 (and 5" — >■ A+p in Example 6), we integrated over the 
time dependence, and took sum over 5° (-B'^) decay events and their CP conju- 
gate. However, if we take difference between the decay events and their CP 
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Imi^C^C;]), (54) 

where Xb = Atub/Tb- If we detect at least one of them, it means that CP is 
violated. 

In the Examples 1-5, we showed some simple processes. We can determine 
the parameters more precisely by using not only these processes but also r"*" — >■ 
tt'^tt'^Dt and other processes. 

In Section [SI we studied the case that only one of two fermions is unstable. 
B ^ Tfj, mode is new physics itself. So if anything, value is very significant 
for understanding it. 

In Section [6l we studied the baryon modes. These modes contain the non- 
perturbative QCD effects to pollute the possible new physics effect. If we sup- 
pose that there is no new physics, we can test the factorization and the pole 
model, for example. On the other hand, if the experimental result highly differs 
from these predictions, we should consider the new physics effect. Recently, 
the lattice gauge theory predicts some B decay processes [25]. We hope that 
the lattice gauge theory predicts precisely the B meson baryonic decays in near 
future. If so, we can search for new physics, precisely. 

We emphasize that, to discover new physics, it is necessary to determine as 
many physical quantities as possible, and compare them to the SM predictions. 
Moreover, it is preferable to be done by the unified form for simplicity, facility, 
and practicality. This paper will help this process. 



conjugate instead of sum, we obtain 

i + A^ i + A^ q 

i + Ag ^ + ^B 1 
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A differential branching ratios for fa and fb de- 
cays 

The differential Branching ratio for the process fa^o, + anything and fi, 
b + anything in fa and ft rest frame, respectively, are [25] : 



dBr{fa -)■ g + anything) ^ 2 
dBr{fii b + anything) 



d^h 



Bn 



2 

7rm?A6 



(55) 



where G"' {ya.b) and G2 iUa.b) are the functions of ya^b = "^Ea^b/rria.b, K,b are 
defined as 



Aa.6 = / dya.bvl.bGi '' {ya.b) 



(56) 



ka h = ka,fc/|k(j_fc| where ka,b are the momentum of the particle a and 6, respec- 
tively. 

We note here that physical vector quantities which we treat in this process 
are only s"''' and ka,f,. The only scalar made by these vector quantities is s^^^-ka^f, 
So, we can explain the differential branching ratio, Eq. (j55p by only two terms 
which are proportional to Gi'^{ya,b) and s°-'^ ■ ha,bG2{ya,b), respectively. 
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